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In this paper we study the effects of f(R) Theories of Gravity on Solar System gravitational 
tests. In particular, starting from an exact solution of the field equation in vacuum, in the Palatini 
formalism, we work out the effects that the modifications to the Newtonian potential would induce on 
the Keplerian orbital elements of the Solar System planets, and compare them with the latest results 
in planetary orbit determination from the EPM2004 ephemerides. It turns out that the longitudes 
of perihelia and the mean longitudes are affected by secular precessions. We obtain the resulting 
best estimate of the parameter k which, being simply related to the scalar curvature, measures the 
non linearity of the gravitational theory. We use our results to constrain the cosmological constant 
and show how f(R) functions can be constrained, in principle. What we obtain suggests that, in 
agreement with other recent papers, the Solar System experiments are not effective to set such 
constraints, if compared to the cosmologically relevant values. 
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I. INTRODUCTION 

In order to make General Relativity (GR) agree with 
the recent observations in Cosmology, which evidence the 
acceleration of the Universe thanks to experimental data 
coming from different tests [2, H, H, 0, [H] , the existence 
of the dark energy is needed. On the other hand, the- 
ories alternative to GR, such as Higher Order Theories 
of Gravity (HOTG) can explain the acceleration of the 
Universe without requiring the existence of dark energy 
[I, 0, H H. The simplest model of HOTG are the so- 
called f(R) theories, where the Lagrangian depends on 
an arbitrary analytic function / of the scalar curvature 
R. As it is well known, f(R) theories can be studied both 
fmjill, [H [H and in the Pala- 
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tini formalism [ll [IS EalHllaTialia S3 • We remind 
that, in the vacuum case, f(R) theories, in the Palatini 
formalism, are equivalent to GR with a (non dynamical) 
cosmological constant. 

Since GR is in sharp agreement with the experimen- 
tal results obtained in the Solar System [z| , every theory 
that aims at agreeing with experimental results at cosmo- 
logical scale, should reproduce GR at the Solar System 
scale. 

So, the problem arises of testing the reliability of these 
theories with Solar System experiments. The dynami- 
cal equivalence between f(R) and scalar-tensor theories 
of gravity (see for instance [231 ] and references therein) 
has been used to set bounds on the analytic form of the 
functions f(R), thanks to post-Newtonian parameters of 
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scalar-tensor gravity [2J] . However, this approach was 
recently criticized [251 ] . since the dynamical equivalence 
has been used beyond its range of validity. Actually, the 
debate is still open, and the problem can be faced both in 
the metric formalism and in the Palatini formalism (see 
[2SH3 and references therein). 

Also the issue of the dynamics of f(R) theories in both 
formalisms in the presence of matter (i.e. inside the 
sources of the gravitational field) is under debate [28ll29j|. 
In particular, a recent paper [30( states that f(R) theo- 
ries where a 1/R term is added to the Einstein-Hilbert 
action do not match Solar System experiments: these 
conclusions have been subsequently generalized to other 
forms of f(R) [3Xl ] . Indeed, these analyses are carried 
out in the metric formalism and lead to the same re- 
sults obtained by Chiba [28|. The difficulties arising in 
the metric formalism with the Solar System experiments, 
which ultimately depend on the matching between the 
solutions inside and outside the matter distribution (i.e. 
the matching between the star interior and the vacuum 
exterior), are not present in the Palatini formalism, as 
it has been recently showed 32]: in the Palatini formal- 
ism of the f(R) theories, the space-time inside a star 
does not affect the space-time outside it, (i.e. the vac- 
uum solution), contrary to what happens in the metric 
formalism, even though in this approach the mass of the 
star and its density have a non standard relation. 

In a previous paper [26| . one of us studied the con- 
sistence of f(R) theories with observational data in the 
Palatini formalism. We aimed at understanding the cor- 
rections to GR arising from specific modifications of the 
Hilbert-Einstein Lagrangian. We found an exact solu- 
tion (corresponding to constant scalar curvature R) to 
the field equations in vacuum and we showed that the 
modifications to standard GR are directly related to the 
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solutions of what we called the structural equation. The 
latter is an algebraic scalar-valued equation, that con- 
trols the solutions of the field equations of the theory. 
These modifications can be suitably interpreted as post- 
Newtonian parameters related to the non linearity of the 
theory. 

The exact solution we found in 26] corresponds to the 
Schwarzschild-de Sitter metric, which has been studied 
in the past [H, [H, [H[ and, also, more recently in con- 
nection with the problem of the cosmological constant 
[36l . [37 . 38] . The relevance of the Schwarzschild-de Sitter 
metric on f(R) theories of gravity has also been discussed 
in [3^]. Here, we explicitly work out the effects of f(R) 
theories on the Keplerian orbital elements of the Solar 
System planets, and compare them with the latest re- 
sults in planetary orbit determination from the EPM2004 
ephemerides, in order to obtain the best estimate for the 
parameter k, which is a measure of the non linearity of 
the theory. We compare this estimate to the values of the 
cosmological constant A and suggest that it can be used 
to set bounds on the parameters of the f(R) functions 
relevant in cosmology. 

The paper is organized as follows: after briefly intro- 
ducing the theoretical framework of f(R) gravity and the 
vacuum exact solution in Sec. HH we investigate the per- 
turbations of the orbital elements in Sec. IIII1 which are 
compared to recent data in Sec. IIVI then, we show how 
our approach may lead to constraints in cosmology in 
Sec. El Conclusions are outlined in Sec. IVII 



II. 



VACUUM EXACT SOLUTION OF f(R) 
GRAVITY FIELD EQUATIONS 



The equations of motion of the f(R) gravity in the 
Palatini formalism can be obtained, by independent vari- 
ations with respect to the metric and the connection, 
from the action [53 



A = A grav + A mat = / [y/gf(R) + 2xL mat (ip, Vrp)} d 4 x, 

(1) 

where R = R(g,T) = g a ^R a0 (T), R^(T) is the Ricci- 
like tensor of any torsionless connection T independent 
from the metric g, which is assumed here to be the phys- 
ical metric. The gravitational part of the Lagrangian 
is represented by any real analytic function f(R) of the 
scalar curvature R. The total Lagrangian contains also 
a first order matter part £ ma t functionally depending on 
yet unspecified matter fields ^ together with their first 
derivatives, equipped with a gravitational coupling con- 
stant x = (see e.g. 0,53). 

According to the Palatini formalism [3, HH, OS Ell EE 
EI 0, El, S El H], from Q we obtain the following 
equations of motion: 



1 



f'(R)R (tw) (T)--f(R)g^ 



rpmat 



(2) 
(3) 



where T/^ at — — -2= S g™ tLt denotes the matter source 

stress-energy tensor and V r means covariant derivative 
with res pec t to the connection T. Actually, it is possible 
to show [43LI441] that the manifold M, which is the model 
of the space-time, can be a posteriori endowed with a 
bi-metric structure (M, g, h) equivalent to the original 
metric-affine structure (M, g,T), where T is assumed to 
be the Levi-Civita connection of h. The two metrics are 
conformally related by 



h^v = f'(R) g^v. 



(4) 



The equation of motion ([2]) can be supplemented by the 
scalar- valued equation obtained by taking the g-trace of 
©, where we set r = trT = g^ v T™ at : 



f'(R)R-2f(R) = X r. 



(•5) 



The algebraic equation §E§ is called the structural 
equation and it controls the solutions of equation (|2|). 

The field equations ([2][3]) and the structural equation 
(|5|) in vacuum become 

[f(R)]R { ^(r)-^[f(R)]g^ = 0, (6) 

<(V9[f'(R)]gn - 0, (7) 

f'(R)R-2f(R) = 0, (8) 

As shown in [26[ (see, in particular, Section 3), the system 
of equations ([5][5]) has the spherical symmetrical solution 



) / 2tti kr^ \ > 
ds 2 = - [ 1 - — + — ) dt 2 

+ r 2 dd 2 + r 2 sin 2 tfefy? 2 , 



dr 2 



(9) 



where m is the mass of the source of the gravitational field 
and k = Cj/4, where R — a is any of the solutions of the 
structural equation l[8|). In doing so, we have obtained a 
solution with constant scalar curvature R. 

In particular, if f(R) = R, i.e. our theory is GR, 
R = is the solution of the structural equation, and © 
reduces to the classical Schwarzschild solution. Again, 
for a thorough discussion we refer to [26| , and references 
therein. 

The metric (0) corres pon ds to the Schwarzschild-de 
Sitter space-time (see [45ll46|). which is exactly a solution 
of standard GR with fixed value cosmological constant. 

From © it is evident that the modifications to the 
standard Schwarzschild solution of vacuum GR are pro- 
portional to the values of the Ricci scalar, owing to the 
proportionality between k and Cj. Indeed, this contribu- 
tion to the gravitational potential should be small enough 
not to contradict the known tests of gravity. In the cases 
of small values of R (which surely occur at Solar Sys- 
tem scale) the Einsteinian limit (i.e. the Schwarzschild 
solution) and the Newtonian limit are recovered, as it is 
evident from ([S]), and it has been proved in [26]. This 
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implies, also, that in order to study the corrections de- 
riving from f(R) theories on the orbital elements in the 
Solar System, the terms proportional to k can be treated 
as a perturbation. 

For further convenience, we remark that, thanks to the 
following change of the radial variable 



r = r I 1 



to kr 

the metric © can be written in isotropic form 



(10) 



ds A = -\l 



2m kr 



dr 



1 + 2m + kr*\ ,^_ 2 + _ %M% + _ 2 s . n2 2) 
r 6 / 

(11) 



up to first order in m and fc. 



III. THE IMPACT OF f(R) GRAVITY ON THE 
PLANETARY ORBITS 



According to what we have seen before, the modifica- 
tions to the solutions of the field equations due to f(R) 



theories are given by a term proportional to the Ricci 
scalar. As a consequence, from (TlTj) we can consider a 
perturbation of the gravitational potential in the form 



AU = Kr 2 , 



(12) 



where k = fc/3. In doing so, we neglect the effect of 
spatial curvature. 

From the potential (fT2"|) we obtain an entirely radial 
acceleration 



-2nr. 



(13) 



Its effect on planetary motions, which is, of course, 
much smaller than usual Newtonian gravity, can straight- 
forwardly be calculated within the usual perturbative 
schemes (see, for instance, [47|). i.e. using the Gauss 
equations, which enable us to study the perturbations 
of the Keplcrian elements, induced by generic perturb- 
ing accelerations. Namely, the Gauss equations for the 
variations of the scmimajor axis a, the eccentricity e, the 
inclination i, the longitude of the ascending node fi, the 
argument of pericentre to and the mean anomaly M. of a 
test particle in the gravitational field of a body to are, in 
general, given by 



da 
It 
de 
dt 
di 
dt 

dn 

~dt 
dio 
~dt 
dM 
~dT 



n\J\ — e 2 



eA r sin ip + A t ^— ^ 



A r sin ip + A t 



COS if - 



e V a) 



i hi 



= A n (-) CQS(W + ip), 



1 



na smi 
— cos i 



Wl — e s 

m 

dt 



A n sin(cj + ip), 



dn y/T^e* 



-A r cos ip + At [ 1 H — | Mil p 



n A r - - VI - e 2 — + cosi— , 

na \a/ \ dt dt 



d.aj 



dn 



(14) 
(15) 
(16) 
(17) 

(18) 
(19) 



in which n = 2ir/P is the mean motion [5£|, P is the test 
particle's orbital period, p is the true anomaly counted 
from the pericentre, p = a(l — e 2 ) is the semilatus rec- 
tum of the Keplerian ellipse, A r , A t , A n are the ra- 
dial, transverse (in-plane components) and normal (out- 
of-plane component) projections of the perturbing accel- 
eration A, respectively, on the frame {r,t, n} comoving 
with the particle. In our case, the perturbing accelera- 
tion p3[) is entirely radial and, consequently, the Gauss 
equations reduce to 



dt 
duj 
~dt 
dM 
~dT 



2c 



da 



dt ?iVl - e 2 



A r sin p, 



na 



de x/T - ? 



dt 
di 

dt=°> 

? = o, 



A r sin<^, 



VT 



nae 
2 



A r cos p, 



n A 

na 



(20) 

(21) 
(22) 
(23) 
(24) 
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TABLE I: Determined extra-precessions Ari7 b s of the longi- 
tudes of perihelia of the inner planets, in arcseconds per cen- 
tury, by using EPM2004 with the default values (3 = 7 = 1, 
= 2 x 10 -7 . The gravitomagnetic force was not included in 
the adopted dynamical force models. Data taken from Table 
3 of H^]. It is important to note that the quoted uncertainties 
are not the mere formal, statistical errors but are realistic in 
the sense that they were obtained from comparison of many 
different solutions with different sets of parameters and obser- 
vations (Pitjeva, private communication 2005a). The corre- 
lations among such determined planetary perihelia rates are 
very low with a maximum of about 20% between Mercury and 
the Earth (Pitjeva, private communication 2005b). 

Mercury Venus Earth Mars 

-0.0036 ± 0.0050 0.53 ± 0.30 -0.0002 ± 0.0004 0.0001 ± 0.0005 



As a result, the inclination and the node are not per- 
turbed by such a perturbing acceleration. 

By evaluating eq. (|13[) on the unperturbed Keplerian 
ellipse 



a(l-e 2 ) 

r = , 

1 + e cos if 

inserting it into (|20j) - (|2"5|) and averaging them with 
dt 1 (1 - e 2 ) 3 / 2 



P 2ir (1 + e cos if) 



7 d<p 



(26) 



(27) 



one gets that only the argument of pericentre and the 
mean anomaly are affected by extra secular precessions 



3 k 



if) = 

at n 



UM\ 
\dt I 



3k (1 



(28) 
(29) 



As a consequence, also the mean longitude A = fi + lu + 
M, which is used for orbits with small inclinations and 
eccentricities, as those of the Solar System planets, is 
affected by an extra secular rate 



I d\\ _ 3k (1 
\dl] ~ ~n 



y/T 



(30) 



Among them, there are certainly the totally unmodelled 
general relativistic gravitomagnetic field, which induces 
the Lense-Thirring planetary precessions, and the solar 
quadrupole mass moment which, instead, was in- 
cluded in EPM2004, but it is currently affected by a 
~ 10% uncertainty. For a discussion of such issues in 
a different context see [50]. Each planetary perihelion 
is, thus, affected by such residual effects, so that it is 
not possible to entirely attribute the determined extra- 
advances to the action of AU of eq. (fT2")) . especially for 
Mercury for which the mismodelled/unmodelled Newto- 
nian and Einstcinian perturbations are stronger. A bet- 
ter way to use Am is to suitably combine them in order 
to make the estimate of k independent, by construction, 
,of the Lense-Thirring an d J® effects. According to the 
approach followed in l5l| . it is possible to construct 
the following combination 



Am 



Mc 



biAm 



• Ear 



b 2 Am 



Ma 



-1.97 x 10 9 s 



(31) 



with bx = —80.7 and b 2 = 217.6. The dimensionless coef- 
ficients 61 and b 2l which are built up with a, e and i of the 
planets adopted, cancel out, by construction, the impact 
of the gravitomagnetic field and of the solar quadrupo- 
lar mass moment on the combination eq. (|31|) . This 
can straightforwardly be checked by combining the per- 
ihelion precessions due to J 2 and to the Lense-Thirring 
force with the coefficients of eq. (|3"Tj) : the result is zero. 
Eq. (jni|) and the values of Table [fl converted in s , 
allow to obtain 



-2.6 x 10" 26 s~ 2 ' 



(32) 



The error on k obtainable from eq. (|3ip can conserva- 
tively be evaluated as f5£j 



8k < 



8(Am 



Mc 



fei|<5(Ari7 tar ) + \b 2 \8{Am 



Ma 



1.9 x 10 9 s 



1 x 10~ 25 s~ 2 . 



(33) 



V. DISCUSSION 



Because of the relations k — 3k, and k — R/A, from 
the estimate (|3"2"]) on k, we get corresponding bounds on 
the allowed values of the scalar curvature R: 



IV. COMPARISON WITH THE LATEST DATA 



The expressions obtained are useful for compari- 
son with the latest data on planetary orbits from the 
EPM2004 ephemerides [48| by exploiting the determined 
extra-advances Am of the longitudes of perihelia m of 
the inner planets (Table 3 of [49j, here partly reproduced 
in Table©. 

The extra-advances Am determined in (49| are af- 
fected, in general, by all the Newtonian and non- 
Newtonian features of motion which have not been ac- 
counted for in the dynamical force models of EPM2004. 



\R\ < 3.12 x 10" 



(34) 



or, restoring units in such a way that R is measured in 
2 we have 



\R\ < 3.47 x 10 



-42 -2 

m . 



(35) 



Since the k parameter of the Schwarzschild-de Sitter 
metric can be interpreted as a cosmological constant A 
(see [1^]), we get for A the following limits 
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\k\ = |A| < 8.68 x l(T 43 m- 2 . (36) 

This result is a much smaller than the one obtained in 
36|. and comparable to the those obtained in [37l | and 
52], but it is still several orders of magnitude greater 

than the current value of the cosmological constant Ao ~ 

10" 52 m- 2 m. 

As we have seen in Sec. [Til the k parameter is simply 
related to the solutions of the structural equation 

f'(R)R-2f(R) = 0, (37) 

namely it is k = Ci/A, where R = Ci are the solutions of 
eq. (I37p . As a consequence, the bounds on the n param- 
eter, in principle, enable us to constrain the functions 
f(R). We may proceed as follows. 

In general, the functions /, beyond the scalar curvature 
R, depend on a set of N constant parameters aj, j = 
1..N, so that we may write / = /(i?, ai, .., aw), and, on 
solving eq. (|3"T)) . we obtain 

R = T(ai, .., ajv). (38) 

Consequently, what we ultimately obtain is a limit on 
the allowed values of the combination T of these param- 
eters. 

In Q the Lagrangian 

f(R)=R-^ (39) 

was introduced, and it was proved that it mimics cos- 
mic acceleration without need for dark energy. However, 
we point out that the Lagrangian (|39[) was found to be 
instable [l(| El, H3|- For discussion on this issue, we re- 
fer to the recent paper by Faraoni [55|. Nonetheless, we 
consider the Lagrangian (|39p as an example in order to 
illustrate our approach, while a more general discussion 
will be carried out elsewhere [561 ] . 

The Lagrangian (|39|) depends on the parameter /i only, 
and eq. (f3"8"]) becomes 

\R\ = V3fJ?, (40) 

and, thanks to ()36|) . we may set a limit on the parameter 
fx: 

fi < 1.41 x 10~ 21 m _1 , (41) 

or, since \x is a parameter whose physical units are those 
of a mass, 

fi < 2.80 x 10" 28 eV. (42) 

This value is remarkably greater than estimate y, ~ 
10~ 33 eV @, needed for f(R) gravity to explain the ac- 
celeration of the Universe without requiring dark matter. 

The same result hold for the slightly different La- 
grangian f(R) = R- l j T + ai? 2 . 

The same approach that we have briefly outlined here, 
can be applied to other functions f(R), and it could be 
useful, in particular, for those for which the structural 
equation can be solved analytically. 



VI. CONCLUSIONS 

In this paper, we have studied the f(R) Higher Or- 
der Theories of Gravity and, in particular, the effects 
of the modifications of the Newtonian potential on the 
Keplerian orbital elements of the Solar System planets. 
Starting from an exact solution of the field equations in 
vacuum in the Palatini formalism, which corresponds to 
the Schwarzschild-de Sitter metric, we have considered 
the perturbations of the orbital elements, by means of 
the Gauss equations. We have showed that the longi- 
tudes of perihelia and the mean longitudes are affected by 
secular precessions and, by comparison with the latest re- 
sults in planetary orbit determination from the EPM2004 
ephemerides, we have obtained the resulting best esti- 
mate of the parameter k which measures the non linearity 
of the gravitational theory. Since the k parameter can be 
interpreted as a cosmological constant A, our approach 
enables us to set the limit |A| < 8.68 x 10~ 43 m~ 2 , which 
is much smaller than the one obtained in [36j and com- 
parable to the those obtained in [37], [HJ , but still too big 
if confronted with the current value of the cosmological 
constant Ao ~ 10 _52 m -2 . Furthermore, we have showed 
that the best estimate of the parameter k, can be used 
to constrain the functions f{R). In particular, our ap- 
proach is suitable for those f(R) for which the structural 
equation (|3"T|) can be solved analytically, thus allowing an 
explicit evaluation of the parameters appearing in them 
from cosmological experiments results. As an example, 
we have considered f{R) = R — y 4 /R, and we have ob- 
tained the limit for the parameter /i : /i < 2.80 x 10 _28 eV, 
much greater than /i ~ 10~ 33 eV, required in order to 
agree with cosmological observations [6j. 

What our results, in agreement with [H, H3, [H, 
suggest is that, in general, Solar System experiments are 
not able to constrain k or the other related parameters 
up to orders of magnitude comparable to the cosmolog- 
ically relevant values. On the other hand, this fact can 
be also interpreted by saying that f(R) theories are vi- 
able on the Solar System scale, since their predictions 
are indistinguishable from GR ones, while they are ap- 
preciably different only on much larger scales, such as the 
cosmological one. 

All we have done, refers to a solution of the field equa- 
tion in vacuum, and corresponds to a space-time metric of 
constant scalar curvature R. Actually, a complete study 
of the comparison between GR and f(R) theories needs 
more general situations, such as, for instance, those re- 
quiring solutions of the field equations within the matter 
distribution: this will be done elsewhere. 
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Let the signature of the 4-dimensional Lorentzian man- 
ifold M be ( — ,+,+,+); furthermore, if not otherwise 
stated, we use units such that G = c = 1. 
For an unperturbed Keplerian ellipse it amounts to n = 
^/Gm/a s . 

Although rather small, a certain amount of correla- 
tion among the determined extra-advances of perihelia 
is present, with a maximum of 20% between the Earth 
and Mars (E.V. Pitjeva, private communication, 2006). 
If a root-sum-square calculation is performed, an upper 



bound k = 8 x 10 



is obtained. 



